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ABSTRACT 



Context. The symbiotic system HM Sagittae consists of a Mira star and a secondary White Dwarf component. The dust content of the system 
was severely affected by the nova outburst in 1975, which is still ongoing. The capabilities of optical interferometry operating in the mid-IR 
allow us to investigate the current geometry of the dust envelope. 

Aims. We test our previous spectro-interferometric study of this system with new interferometric configurations, increasing the uv coverage 
and allowing us to ascertain the appearance of the source between 8 and 13 /im. 

Methods. We used the MIDI instrument of the VLTI with the unit telescopes (UTs) and auxiliary telescopes (ATs) providing baselines oriented 
from PA=42° to 127°. The data are interpreted by means of an elliptical Gaussian model and the spherical radiative transfer code DUSTY. 
Results. We demonstrate that the data can be reproduced well by an optically thick dust shell of amorphous silicate, typical of those encountered 
around Mira stars, whose measured dimension increases from 8 to 13 yum. We confirm that the envelope is more extended in a direction 
perpendicular to the binary axis. The level of elongation increases with wavelength in contrast to our claim in a previous study. 
Conclusions. The wider uv coverage allows us to deepen our previous investigations of the close circumstellar structure of this object. 
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1. Introduction 

HM Sge is a D-type (dust dominated I R) symbiotic system 
that erupted as a symbiotic nova in 1975 ( Dokuchaeval 1 1 976h . 



evolving from a 17 to a 11 magnitude star with a rich 
emission-line spectrum. The cool component is a Mira star and 
the hot component is a White Dwarf (WD), which had escaped 
detection being enshrouded in the dense Mira envelope before 
this dramatic event. 
Previous works 



by iSchild et alj (1200 Ih and 



Bogdanov & Taranoval ( 2001 ) were published on HMSge 
the same year, evaluating the parameters of the dust in the 
system by fitting the ISO/SWS data using the DUSTY radiative 



* Based on observations made with the Very Large Telescope 
Interferometer at Paranal Observatory under programs 075.D-0484, 
077.D-0216, and 079.D-0213 

** Reduced visibilities and differential phases are avail- 
able in electronic form at the CDS via anonymous 
ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http : //cdsweb . u- strasbg . f r/cg i -bir\/qcat? J/A+A/?/? I 



transfer code (llvezic et al.L Il999h . Both groups used single 
and double component models but did not reach identical 



conclusions. Schild et al 



(120011) favored a 2-shell model, 
which a ccounted for the effect of the W hite Dwarf on the dust, 
whereas Bogdanov & Taranoval ( 2001 ) proposed a single-shell 
model. Since the extensions of the dusty shell differ from one 
model to another, mid-infrared interferometry data enabled 



Sacuto et al.l (120071) to determine the most suitable of the 
proposed models. The authors showed that the morphology 
of the dust was accounted for well by the single shell model 
of ISchild et al] (120011) . which included a compact and opti- 
cally thick circumstellar environment, typical of a Mira star 
However, because of its high optical depth, the model is 
inconsistent with the prese nce of non-absorbe d emission of 
the Mira in the near-IR, and lSacuto et al.l (l2007h indicated that 
this issue could be solved by assuming a more complex spatial 
distribution of dust. We attempt to confirm the morphology of 
the dusty envelope using observations acquired with five new 
interferometric baselines secured in 2007. 
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2. Observations 

The Very Large Telescope Interferometer (VLTI) of ESO's 
Paranal Observatory was us ed with MIDI, the MID-infrared 



Interferometric recombiner jLeinert et all 120031) . MIDI com 



bines the light of two telescopes and provides spectrally 
resolved visibilities in the N band atmospheric window. 

The first observations of HM Sge were conducted in 2005 
and 2006 with the VLT unit telescopes (UTs) UT2, UTS, and 
UT4, providing projected baselines and projected angles, in 
the range of 32-59 meters and 42° to 105° respectively (see the 
first two pa rts of Table [p. An interpretation of these data is 
described in Sacuto et al. (12007). 



The acquisition images enabled us to confirm that the 
mid-IR source was unresolved by the data from the individual 
UT. This implied that most of the mid-infrared flux originated 
in the inner 300 mas of the source, corresponding to the Airy 
pattern of the UTs and defining the Field Of View (FOV) of 
the interferometric observations. This also means that the flux 
detected within the 1.1" FOV of the VLTI auxiliary telescopes 
(ATs), originates in this inner region. Therefore, the visibilities 
extracted from the ATs and the UTs are fully consistent with 
each other. 

Complementary observations were completed in 2007 with 
the ATs AO, DO, Gl, and HO, providing projected baselines 
and projected angles, in the range of 64-89 meters and 72° 
to 127° respectively (see the third part of Table [TJ. With 
longer projected baselines, these configurations allow us to 
probe smaller circumstellar regions. The new projected angles 
provide additional constraints about the morphology of the 
dust shell. 

Table[T]presents the journal of interferometric observations. 
The calibrators, HD188512 (G8IV, diam=1.98+0.02mas), 
HD187642 (A7V, diam=3.22+0.01 mas), HD187076 
(M2II, diam=8.05+0.13mas), HD 1676 18 (M3.5III, 
diam=11.33+0.04mas), and HD206778 (K2Ib, 
diam=8.38+0.09masU were observed immediately be- 
fore or after each science target observation. The third part 
of the table gathers the additional observations of the object 
obtained in April, May, and July 2007. 

Chopped acquisition images were recorded (f=2Hz, 2000 
frames, 4 ms per frame) to ensure the accurate acquisition 
of the target. The acquisition filter was an N-band filter 
Photometry was obtained before and after the interfero- 
metric observations with the HIGH-SENS mode of MIDI, 
using the prism that provides a spectral dispersion of about 
30. The dat a reduction software package^ MIA and EWS 
(|jaffeLl2004 we re used to prepare the spe ctra and visibilities 
dChesneau et all 120051: iRatzka et alll2007l) . MIA is based on a 
power spectrum analysis and uses a fast Fourier transformation 
(EFT) to calculate the Fourier amplitude of the fringe packets. 



' MIDI calibrator database 
^ http : //www . mpia-hd . mpg . de/MIDISOFT/ , 
|http : //www . strw . leidenuniv . nl/~nevec/MIDl71 



Table 1. Journal of all available MIDI observations. The cali- 
brators used to calibrate the visibilities are given below the sci- 
ence target. The phase of the Mira (<^Mira) during the observa- 
tions is indicated. The configuration used for the observations 
is given. The length and position angle of the projected baseline 
are also indicated. 



Star 



UT date & Time 



Basc[ra] PA[deg] 



HMSgc 
HD188512 

HMSgc 
HD188512 

HMSgc 
HD188512 

HMSgc 
HD188512 



2005-07-23 06:48:42 
2005-07-23 07:12:10 
2005-07-2402:41:11 
2005-07-24 02:58:55 
2005-07-2403:42:18 
2005-07-24 04:03:14 
2005-07-2406:14:49 
2005-07-24 06:37:58 



46.5 



44 



HMSgc 2006-05-17 09:04:02 0.31 

HD187642 2006-05-1709:27:14 

HMSgc 2006-06-11 09:05:09 0.36 

HD187642 2006-06-11 08:42:52 



HMSgc 
HD187076 

HMSgc 
HD167618 

HMSgc 
HD187076 

HMSgc 
HD206778 

HMSgc 
HD206778 



2007-04-04 09:39:22 
2007-04-04 09:06:17 
2007-05-08 08:52:55 
2007-05-08 10:16:28 
2007-05-10 08:34:58 
2007-05-1008:55:17 
2007-05-16 09:48:33 
2007-05-1609:20:11 
2007-07-30 05:54:54 
2007-07-30 05:23:12 



0.93 



while EWS uses a shift-and-add algorithm in the complex 
plane, averaging appropriately modified individual exposures 
(dispersed channeled spectra) to obtain the complex visibility. 
An extensive study of the fluctuations of the instrumental 
transfer function was carried out in all periods. The rms scatter 
of these calibration measurements was ~0.05 between 8 jum 
and 13 /im, corresponding to a relative error of 10%. 

All visibility and differential phase data as well as all the 
characteristics of the observations are available from the CDS 
(Centre de Donnees astronomiques de Strasbourg); all data 
products are stored in the FITS-based, optical-int erferometry 



data-e xchange-format (OI-FITS) described by iPauls et al 
(l2005h . 



In Fig. [U the MIDI spectra for each period (2005-2006 



2007) is compared with the ISO/SWS spectra dSloan et al 



2003h acquired at two different phases of the Mira. The level 
of flux received by MIDI is equivalent to that of ISO, implying 
that most of the N flux is concentrated within the 300 mas beam 
of a 8 m telescope. It also means that the level of emission has 
probably not changed significantly since the ISO observations 
(1996-1997), keeping in mind a ty pical error for the abs olute 
SWS flux calibration of 15 to 30% dSchaeidt etallll996h . For 
physical and technical reasons, the flux error bars of the pho- 
tometry from the UTs are below the 10% level, whereas the 
accuracy of the photometry extracted from the ATs is much 
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Table 2. Parameters for HMSge deduced from the best-fit 
DUSTY model to the spectrophotometric and interferometric 
data. 



9 10 

wavelength {}irn) 



Fig. 1. MIDI and ISO/SWS flux of the source in the mid- 
infrared. Dotted line with -10% error bars is the 2005 MIDI 
flux from the UTs (0~O.75); dashed line with ~10% error bars 
is the 2006 MIDI flux from the UTs (0-0.30); dashed dot dot 
dot line with -25% error bars is the 2007 MIDI flux from the 
ATs (0-0.00). The two soHd lines ai-e the ISO/SWS spectra 
(upper one: 1996, 0=0.65; bottom one: 1997, 0=0.08) 



worse with an error bar of 25% (see Fig. [TJ. The error levels 
are estimated from the standard deviation of the flux-calibrated 
spectra determined from each night and each telescope, i.e typ- 
ically 4 to 8 spectra. As a consequence, it is difficult to dif- 
ferentiate between photometric variations in one phase from 
another. The number of measurements acquired in each period 
is also limited. Therefore, we decided in the following to con- 
tinue the spectro-interferometric study of the star neglecting its 
variability. 

3. Radiative transfer modeling 

The spectro-interferometric observations are interpreted in 
terms of spherical models generated by the DUSTY package 
dlvezic et al. 11999 '). We first fit the spectrophotometric data 
of the star, then we determine the corresponding synthetic 
spectrally-dispersed visibility profiles throughout the N band 
(8-13;um) and compare them with the MIDI visibilities for 
each baseline. 



3.1. The compact single-dust-shell model 

T he best-fit models ob tained in our previous work (see Sect. 5 
of lSacuto et al.l (l2007n ) involved only a single shell component 
around the Mira. The modeling w as based on the sim ultaneous 
fit of the ISO/SWS spectrometric dSloan et al.[l2003h and inter- 
ferometric data for periods 2005 and 2006 (see the first two part 
of TablelTJ. Different models were obtained: one involving pure 
amorphous silicates, the other including carbon material in the 
dusty layers. We tested these models upon the new MIDI in- 
terferometric measurements (see the third part of Table [TJ and 
also some near-infrared photometric data not used in the first 
study. 

Near-infrared photometric data provide an additional constraint 



Parameter 



Value 



Effective temperature (K) 3000 

Central star diameter (mas) f .6 

Luminosity (Lg) 5000 

Distance (kpc) f.5 

Shell inner radius: e,„ (mas) 3.6 

Inner boundary temperature (K) 1600 

Grain chemical composition 100% Sil-Ow' 

Density power law coefficient 2.2 

Grain size distribution MRN^ 

Geometrical thickness (Bi„) 1000 

Visual optical depth 25.2 

10 //m optical depth 2.15 
'Sil-Ow stands for 'warm' silicates of Ossenkopf et al. (199^. 
'grain size distribution as described bv .Mathis et al.. C1977.V 



of the radiation originating in the Mira star, whereas the inter- 
ferometric data provide the necessary high angular resolution 
to probe the internal layers and to derive information about the 
morphology of the structure. 

A least square fitting minimization was performed taking into 
account spectrophotometric and interferometric data. This min- 
imization was completed for a large grid of parameter values 
by varying the effective temperature of the star from 2800 to 
3200 K (ATeff=100 K), the inner boundary temperature from 
1400 to 1800 K (ATin=100 K), the density power-law coef- 
ficient from 1.8 to 2.4 (Ap=0.1), and the lOyum optical depth 
from 1.5 to 2.5 (ATio^m=0.05). 

Figure|2]shows the best-fit model allowing a simultaneous fit of 
spectrophotometric and interferometric data of the star Table |2] 
presents the parameters of HM Sge deduced from this fit. 

The results do not differ significantly from the parame- 
ters of the best-fit model deduced in the previous work of 
Sacuto et al.l (l2007h . and reproduce well the near-infrared 



photometric data and the new interferometric data of the star. 
However, the new data enable us to exclude the model includ- 
ing carbon dust proposed to improve the fit. This assumption, 
involving a large fraction of carbon, is not physically justified 
given the spectral type of the dust-forming source, an M-type 
Mira star In the absence of any clear signature of other dust 
species and given the spatial complexity of the source, it is 
extremely difficult to use the dust composition as a parameter 
for improving the fit. 

The interpretation by means of a sin gle dust shell sur- 
rounding the Mira star was first rejected bv iSchild et al.l (1200 lb 
because of its high opacity, which was inconsistent with the 
detection of the Mira photometric cycle in the near-infrared. 
However, interferometric data imply that a simple dust shell is 
the most suitable solution. 

With the help of these new observations, the increase of the uv 
coverage allows us now to investigate the evidence of departure 
of the circumstellar environment from spherical symmetry. 
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Fig. 2. Top-left: Best-fit model of the single dust shell spectrum (thin line) supperimposed on the ATs MIDI flux (error bars), 
the 2 ISO/SWS spectra (upper thick line: 1996, 0=0.65, and bottom thick line: 1997 , 0=0.08), and the ph otometric data (open 
triangles) with (U/B/V/R/J/H/K/IJM/N) data from lYudin et all |l994) and / band fromi chochol et alj (|2004 . Top-right: close-up 
view of the best-fit model to the dust feature. Bottom: the corresponding model visibility (solid line) superimposed on the MIDI 
visibilities (error bars) for the eleven projected baselines. 



4. Morphological interpretation 

The phase shifts or differential phases is information extracted 
from the difference between the phase at a given wavelength 
and the mean phase determined in the full N-band region. It 
provides partial information (depending on the uv coverage) 
about the asymmetry of the source, in addition to the infor- 
mation contained in the visibility amplitudes. As an example, 
strong phase signals were detected in sources harboring dusty 



disks dDeroo et all 120071: lOhnaka et all 120081) . 



Previous results of differential phases obtained bv lSacuto et al 



(120071) indicated no obvious signature of any asymmetry 
in the close environment of the star. Due to the low rms 
(< +5°) of the calibrated differential phases averaged over all 
the projected baselines, new interferometric configurations 
confirm the global centrosymmetry of the object in the vicinity 
of the source (see Fig. [3]). This allows us to justify in the fol- 
lowing the utilization of a centrosymmetric geometrical model. 

The geometrical constrai nts obtained from the data pre- 
sented in ISacuto et al.l (l2007h show that the circumstellar en- 
vironment of the star is elongated perpendicularly to the binary 
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Wavelength (/xm) 

Fig. 3. Calibrated differential phase of the object evaluated 
from all the determined ATs projected baselines (dotted lines). 
Error bars correspond to the rms of the differential phase aver- 
aged over all the projected baselines. The bold solid line is the 
central value of these error bars. 



axis (PA=130° NE; 'Eyres et al."200l'), in the direction of the 
Raman line polarization ( Schmid et aL 2QQQ.). This informa- 
tion was obtained by comparing the half width at half maxi- 
mum of the equivalent one-dimensional Gaussian distribution 
from various projected baseline position angles. With the new 
baselines, we can perform a more robust analysis using a two- 
dimensional fitting in the Fourier space by means of an ellipti- 
cal Gaussian flux distribution, appropriate for studying devia- 
tions from spherical symmetry. 

4.1. Elliptical Gaussian model 

The use of purely geometrical models can help to constrain the 
object morphology. Because the circumstellar envelope is opti- 
cally thick (To.55^rn = 25.2; see Table|2]), a Gaussian brightness 
distribution centered on the star is a natural approximation of 
the flux distribution in the plane of the sky. 

The visib ility of th e elliptical Gaussian model can be ex- 
pressed by CTvcner et aL .2004) 



V(u', v') - exp 



;rXj 



Aln2 



(1) 



where flmj is the angular size FWHM of the major axis, r is 
the axial ratio of the minor to the major axis, and (m',v') are 
the modified spatial frequencies. These frequencies can be ex- 
pressed in terms of the spatial frequencies (u,v) by using a co- 
ordinate transformation 



u — u cos (p — V sm ^ 
v' - usmd) + vcos ( 



(2) 
(3) 
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Fig. 4. Spectrally-dispersed uv coverage of the 2005-2006- 
2007 observations from 8 to 12.5 fim. 

respectively, and A is the observed wavelength. Finally, is 
the position angle (measured east from north) of the major axis. 

By completing a non-linear least squares fitting, based on 
the Levenberg-Marquardt method, we obtained the best-fit el- 
liptical Gaussian model to the data with a reduced of 0.27. 
Figure |4] shows the spectrally-dispersed uv coverage of all the 
observations (see TablefT]). Figure|5]provides the brightness dis- 
tribution of the best-fit elliptical Gaussian model of HM Sge 
between 8 and 12.5 yum. Figure |6] shows the resulting FWHM 
of the major axis, and the ratio of the minor to the major axis 
for the best-fit elliptical Gaussian model as a function of wave- 
length. 

The best-fit model is obtained for an extension of the struc- 
ture perpendicular to the binary axis located at PA=130° NE 



where u - -ycos^p and v - -y sin^p, where Bp and Op are 
the projected baseline and the projected angle on the sky 



dEvres et all 1200 lb as presented in our previous work. The 
length of the major ax is is also approximatively the same as that 
published in lSacuto et al. (2007). However, the level of elonga- 
tion increases from 8 to 12.5 fim, in contrast to results found 
before. This difference can be explained by the wider uv cov- 
erage allowing the geometry of the structure to be more tightly 
constrained. 

5. Conclusions 

The wider uv coverage data has allowed us to confirm that the 
dust envelope around the Mira is described well by an opti- 
cally thick dust shell that consists predominantly of amorphous 
silicates. We also confirmed the longer extension of the circum- 
stellar envelope in a direction perpendicular to the binary axis, 
and the increase in the size of the structure between 8 and 12.5 
jum. It has also been demonstrated that the amount of elonga- 
tion is higher at 12.5 //m and steadily decreases toward 8 fim, 
but the ratio between the two axis remains limited at most to 
0.85, implying only a slight departure from spherical symme- 
try. This behavior can naturally be interpreted. The layers emit- 
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Fig. 5. Brightness distribution of the best-fit elliptical Gaussian model of HM Sge from 8 to 12.5 |/m (north is up and east is left). 
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Fig. 6. Le/f: full width at half maximum of the major axis for the best-fit elliptical Gaussian model as a function of wavelength. 
Right: ratio of the minor to the major axis for the best-fit elliptical Gaussian model as a function of wavelength. 



ting most of the 8 jum flux, located closer to the photosphere, 
are the leas t affected by the White Dwarf, situated far (60 AU; 
Eyres et al.ll2001,) from the Mira. At 13 fim, the interferometer 
probes farther regions that are more affected by the energetic 
but diluted wind originating in the erupting White Dwarf. We 
also note that the new baselines, recorded in 2007, range from 
60 to 89m whereas the previous ones ranged from 32 to 59m 
and probed outer regions of the dusty envelope, which were 
potentially more affected by the White Dwarf. Therefore, we 
expect the amount of elongation to be higher, as reported in 
Sacuto et al. 1(2007'), although the poor uv coverage of the 2005- 
2006 observations does not allow us to estimate it accurately. 

These MIDI/VLTI interferometric observations probe a 
compact region that includes the giant star and its vicinity to 
a radius of about 30 mas, set by the smallest baselines (~32m). 
There is a clear requirement for near and mid-IR high dynam- 
ics, imaging at the diffraction limit of a 8-lOm class telescopes, 
to probe the inner nebula (0.8" wide) alrea dy studied with 



the H ST jEvres et all I2OO1I) and in the radio (iRichards et al. 
1999h . 



studying such a complex source. It is also uncertain whether 
the three telescopes recombiner AMBER, operating in the near- 
IR, can provide a substantial increase in information about the 
morphology of the system, because the (spherical) Mira and 
hot dust in the inner wind still dominate. A significant im- 
provement in the quality of image reconstruction must await 
the advent of the second generation instrument MATISSE 



dLopez etal.[l2006h . 



We are close to the limits of the imaging capabilities of 
the VLTI used with the two telescopes recombiner MIDI for 
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